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PREFACE 

The  work  reported  herein  was  performed  by  the  Arnold  Engineering  Development  Center 
(AEDC),  Air  Force  Materiel  Command  (AFMC)  under  program  element  65807F.  The  results  were 
obtained  by  Calspan  Corporation,  AEDC/Operations,  support  contractor  for  aerodynamic  testing 
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37389-4300,  under  AEDC  Project  No.  0107.  The  AEDC  Project  Managers  were  Maj.  M.  Maclin 
and  Maj.  Mark  Briski.  The  manuscript  was  submitted  for  publication  on  July  6, 1994. 
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1.0  INTRODUCTION 

Exhaust  gases  from  rocket  motors  and  outgassing  products  of  spacecraft  materials  can 
migrate,  in  a  space  environment,  to  other  surfaces  of  the  spacecraft  and  affect  their  functions 
(Ref.l).  Cryogenically  cooled  optical  systems  are  particularly  vulnerable  to  these  contaminants  as 
are  thermal  control  surfaces  and  solar  cells. 

Recent  AEDC  research  has  concentrated  on  the  degradation  of  optical  surfaces  by  condensed 
outgassing  products  of  spacecraft  materials  (Refs.  2, 3,  and  4).  A  prominent  feature  in  the  thin-film 
infrared  transmittance  spectra  of  many  of  these  contaminants  is  a  strong  absorption  band  centered 
at  about  33(X)  cm’’  (Ref.  4).  Other  deleterious  absorptions  frequently  occur  at  wave  number  values 
around  KKX)  cm''  and  less. 

One  of  the  goals  of  the  contamination  program  at  AEDC  is  to  predict  the  optical  effects  of 
contaimnant  layers  on  optical  systems  from  laboratory  measurements  of  the  optical  indices,  n  and 
k,  of  contaminant  materials.  Because  the  cloud  of  contaminant  molecules,  which  eventually  form 
a  condensate,  contains  a  mixture  of  molecular  species,  researchers  are  seeking  ways  to  predict  the 
optica]  properties  of  films  composed  of  constituents  using  the  known  properties  of  each  constituent. 

In  this  report,  the  authors  examine  the  physics  which  describes  the  optical  properties  of  thin 
films  composed  of  mixtures  of  constituents.  We  make  use  of  the  Loientz-Lorenz  relationship 
between  the  optical  indices  of  a  material  and  its  polarizability  to  obtain  the  optical  constants  of  a 
mixture  of  materials.  This  description  is  complicated  by  a  lack  of  detailed  knowledge  of  the  inter- 
molecular  interactions  for  the  molecules  in  a  given  constituent  and  in  a  mixture  of  constituents.  An 
outstanding  problem,  for  example,  is  to  predict  accurately  the  degree  of  hydrogen-bonding  present 
in  any  material  (Ref.  5). 

It  is  also  possible  that,  because  different  materials  have  different  outgassing  rates,  the 
contaminant  layer  is,  actually,  a  “stack”  of  several  thin  films,  each  of  which  has  a  different  compo¬ 
sition.  Interference  effects  within  each  film  may  or  may  not  be  present. 

The  authors  present  several  examples  of  the  analysis  methods  for  thin  films  composed  of 
mixtures  using  AEDC  data  on  20K  thin  films  of  simple  molecules  and  more  recent  AEDC  data  of 
80K  films  condensed  from  the  outgassing  products  of  selected  spacecraft  materials. 


7 


AEDC-TR-94^ 


2.0  LORENTZ-LORENZ  ANALYSIS  METHODS 

2.1  THE  LORENTZ-LORENZ  RELATION 

The  LorentZ'Lorenz  relation  connects  the  complex  optical  index  (n  =  n  -t-  ik)  of  a  material 
with  the  number  density  (N)  of  its  constituent  molecules,  its  mean  polarizability  per  molectile  (a), 
and  a  parameter,  b,  which  relates  the  “local”  electric  field,  E/oc  ,  actually  “felt”  by  a  molecule,  with 
the  electric  field  measured  in  a  cavity  of  the  material,  E .  For  a  linear,  isotropic,  homogeneous  (lih) 
dielectric  material,  the  classical  Lorentz-Lorenz  relation  is: 

Naa  ss  (n2-l)/(n2  +  b),  (1) 

where 

b=(Eoa)“^-l.  (2) 

The  constant  is  the  electric  permittivity  of  free  space,  a  is  another  “local”  field  parameter 
defined  by 

hoc  =  ^  +  a?,  (3) 

and  P  is  the  polarization  density  vector  of  the  dielectric.  If  the  lih  dielectric  contains  only  nonpolar 
molecules, 

a  =  (3  Eq  )  and  b  =  2.  (4) 

(See  Appendix  A  for  a  cursory  derivation  of  the  Lorentz-Lorenz  relation.) 

2.2  CALCULATION  OF  OPTICAL  CONSTANTS  FOR  THE  SAME  MATERIAL  AT 
DIFFERENT  DENSITIES 

One  can  use  the  relation  in  any  of  several  ways.  For  instance,  for  a  single  material  the  Lorentz- 
Lorenz  relation  between  its  optical  index  and  its  (mass)  density,  p  =  mN,  where  m  is  the 
(average)  mass  of  a  molecule,  is, 

p  '  (n^  -  1  )  /  (n^  +  b)  =  aa/m,  a  constant.  (5) 

This  constant  value  (for  a  given  material)  is  often  dubbed  the  “Lorentz-Lorenz  constant”  and 
provides  a  way  to  calculate  the  optical  constants  of  a  material  at  a  density  value  different  from  the 
one  at  which  they  were  originally  measured. 
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23  CALCULATION  OF  THE  OPTICAL  CONSTANTS  OF  A  MIXTURE 

One  can  also  use  the  Loientz-Lorenz  relation  to  calculate  the  optical  constants  of  a  mixture 
of  dielectric  materials  from  the  optical  constants  measured  for  pure  samples  of  each  constituent  ma¬ 
terial.  The  key  property  of  the  Lorentz-Lorenz  relation  that  makes  this  possible  is  that  the  index,n^, 
measured  for  a  pure  sample  of  the  jth  constituent  material,  is  related  to  the  (average)  polarizability 
of  the  molecules  contained  in  the  pure  jth  material  as  follows: 

“j  =  I"  j  -  ^  ‘'j  J*  (6) 

where  Npj  is  the  number  density  of  the  jth  constituent  in  the  pure  sample,  and  aj  and  bj  are  the 
propriate  local  field  parameters  for  the  pure  sample.  Furthermore,  the  (average)  molecular  polariz¬ 
ability  of  the  mixture,  (X,  is  the  weighted  mean  of  the  polarizabilities  of  each  constituent  of  the 
mixture,  Oj;  that  is: 

a  =  lNjaj/N  (7) 

where  is  the  number  density  of  the  jth  constituent  in  the  mixture,  and  N  is  the  total  number  den¬ 
sity  of  the  molecules  in  the  mixture  (N  =  E  Nj).  When  Eqs.  (1 ),  (6),  and  (7)  are  combined,  the  com¬ 
plex  index  of  the  mixture  turns  out  to  be 


n  =  [(l+bA)/(l-A)]^  (8) 

and  A  is  a  complex  number  related  to  the  average  molecular  mass  of  the  mixture  (m  =  £  N^m/N), 
its  local  field  parameters  (a  and  b),  its  density  (p ),  and  the  densities  of  the  pure  materials  (p  pj), 

A  =  [(p/l+b)]Z  [mj/m(l  +  bj)/ppj][ii/-lJ[njVbJ.  (9) 

(See  Appendix  A.) 

Equations  (8)  and  (9)  form  the  heart  of  the  algorithm  of  program  NMIX  which  predicts  the 
optical  constants  of  a  film  material  which  is  a  mixture  of  materials  whose  optical  constants  in  their 
pure  states  are  known. 

3.0  OPTICAL  ANALYSIS  METHODS  FOR  MULTILAYERED  DEPOSITS 

If  the  deposit  on  a  surface  consists  of  a  "stack"  of  uniform  layers,  each  of  which  has  a  distinct 
material  composition,  the  analytical  model  contained  in  AEDC  programs  such  as  TRNLIN  and 
CALCRT  needs  to  be  revised  to  treat  the  deposit’s  optical  properties.  The  analytical  model  of  a 
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multilayered  system  is  often  given  in  terms  of  matrix  mathematics  because  the  effect  of  each  layer 
is  contained  in  its  own  associated  '‘transfer”  matrix.  The  optical  effect  of  the  stack  is  found  by  com¬ 
puting,  in  the  proper  order,  the  product  of  the  matrices  for  each  layer.  It  is,  therefore,  relatively  easy 
to  add  the  effects  of  additional  layers  to  a  given  stack  by  multiplying  appropriate  additional  matrices 
to  the  matrix  of  the  original  stack. 

The  transfer  matrix  of  a  film  which  allows  interference  effects  (i.e.,  allows  phase  coherence 
of  the  beam)  relates  amplitudes  of  the  components  of  the  beam  because  amplitudes  contain  phase 
information.  The  transfer  matrix  of  film  in  which  phase  coherence  does  not  occur  relates  beam 
intensity  components  and  is  quite  different  from  the  matrix  of  a  film  in  which  phase  coherence 
does  occur. 

The  transfer  matrix  of  a  film  is  often  given  as  a  product  of  matrices,  one  or  two  of  which 
represent  the  changes  in  the  beam  at  its  interfaces  and  the  other  gives  the  effects  on  the  beam  as  it 
travels  through  the  film. 

3.1  OPTICAL  MODEL  IN  THE  CASE  OF  PHASE  COHERENCE  IN  THE  HLMS 

In  the  model  we  present  here,  based  on  the  article  by  Potter  (Ref.  6)  and  briefly  discussed  in 
Appendix  B,  the  beam  angle  of  incidence  can  be  from  0  to  90  deg.  This  means  that  each  electro¬ 
magnetic  wave  in  the  beam  has  an  electric  field  with  two  components:  the  electric  field  component 
perpendicular  (TE)  to  the  plane  of  incidence  (defined  by  the  incident  beam  and  the  normal  to  the 
incident  surface)  and  the  component  parallel  (TM)  to  the  plane  of  incidence.  Because  any  electric 
field  can  be  considered  a  vector  sum  of  these  two  components,  we  can  study  the  case  of  a  wave  with 
just  a  perpendicular  electric  component  (the  TE  mode,  or  TEM)  separately  fiom  the  case  of  a  wave 
with  just  its  parallel  electric  component  (TM  mode,  or  TMM). 

3.1.1  TE  Mode 

If  we  choose  a  coordinate  system  where  the  z  axis  lies  along  the  normal  of  the  deposit  (and 
substrate)  surfaces,  the  x  axis  is  perpendicular  to  the  plane  of  incidence,  and,  thus,  the  y  axis  is 
parallel  to  the  plane  of  incidence,  we  can  separate  the  spatial  coordinate  and  time  t  dependencies  of 
the  electric  and  magnetic  vectors  of  each  electromagnetic  wave.  In  the  TE  (perpendicular)  mode, 
the  electric  field  in  a  given  medium  of  index  n  lies  along  the  x  axis  and  has  the  form  of  a  product 
of  functions  of  a  separate  variable, 

E,(y,z,t)  =U(z)X(y,  t),  (10) 


where 
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Uz  =  E^exp  [i2>c'unzcos<|>] ,  (11) 

and 

X(y,  t)  =  exp  [i(27iunysin<ti-a>0]  (12) 

for  an  electric  field  of  amplitude  E^,  vacuum  wave  number  \) ,  angle  from  the  surface  normal,  and 
angnlar  frequency  ci) .  Similar  expressions  occur  for  the  y  and  z  components  of  the  magnetic  field. 

3.1.2  TMMode 

In  the  TM  (parallel)  mode,  the  magnetic  field  is  in  the  x  direction  and  has  the  form, 

Hx(y,z,  t)  =  U(z)X(y,  t),  (13) 

where 

U  (z)  =  Hgexp  [i2i[vnzcos^j  (14) 

for  a  magnetic  field  of  amplitude  Ho,  and  X  is  defined  in  Eq.  (12).  In  this  report  we  follow  the 
convention  that  certain  quantities,  like  U,  have  different  definitions  in  each  mode.  Analogous  to 
the  TE  mode  equations,  there  are  expressions  similar  to  Eq.  (13)  for  the  y  and  z  components  of  the 
electric  field. 

3.U  Phase  Coherent  IVmisfer  Matrix  at  an  Interface  of  Two  Madiii 

Because  the  electric  and  magnetic  fields  of  an  electromagnetic  wave  depend  on  one  other  via 
the  Maxwell’s  Equations,  we  are  able  to  study  the  propagation  of  a  wave  through  a  multilayered 
stack  (in  the  z  direction)  by  noting  what  happens  just  Co  the  value  of  U  in  either  mode.  Furthermore, 
we  can  distinguish  the  portion  of  the  fields  which  propagates  in  the  +z  direction  from  that  part  going 
in  the  -z  direction. 

Suppose  that  U  in  the  Jth  medium  has,  at  the  surface  the  beam  exits,  a  component 
associated  with  the  wave  going  in  the  +z  direction  and  a  component  U^";  let  Uo  J  and  repre¬ 
sent  the  corresponding  quantities  in  the  mth  medium  at  the  surface  the  beam  enters  it.  At  an  inter¬ 
face  between  the  jth  and  mth  media  (Fig.  1),  these  quantities  are  related  by  a  matrix  equation. 
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where  the  amplitude  transfer  matrix  Sj„  is 


Sjm  =  — 

jm  j 

> 


1  r 


jm 

*J 


The  Fresnel  transmission  and  reflection  coe^icients,  and  rj„,  respectively,  are  given  by 

i,.  =  =  i+v  06) 

The  quantity  depends  on  properties  of  the  jth  material  and  on  the  mode  condition  of  the  wave 
(Appendix  B). 


3.1^4  Phase  Coherent  IVansfer  Matrix  Through  a  Uniform  Medium 


The  transfer  matrix  L  for  a  beam  going  through  a  uniform  medium  connects  the  values  of 
both  compqnents  of  U  where  the  beam  enters  the  medium,  and  Ug  with  their  values  a  distance 
d  ficom  that  point  in  the  +z  direction,  and  U  (Fig.  2).  For  a  given  medium,  the  relationships 
among  these  quantities  are 


=  L 

U" 

_u' 

where  the  amplitude  transfer  matrix  is 


L  = 


exp  (-i6)  0 

0  exp  (i5) 
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The  quantity  S  depends  on  the  optical  index  n  and  the  incident  angle  of  the  beam  on  the  first  surface 
of  the  stack  (Appendix  B). 

Note  that  the  quantities  n,  a ,  A,  U ,  if ,  U  ,  S,  t,  r,  ^ ,  L,  and  S  all  contain  phase  information  and 
are,  consequentiy,  complex  numbers. 
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3.1.5  The  Transmittance  and  Reflectance  of  a  Beam  Incident  on  Thin  Films  Deposited  on  a 

Substrate  in  the  Case  of  Phase  Coherency  Within  the  Film  Layers 

A  stack  configuration  often  encountered  experimentally  is  that  of  a  series  of  uniform  films 
deposited  on  one  or  both  sides  of  a  substrate  (Fig.  3).  Using  the  transfer  matrices  in  Eqs.  (15)  and 
(17),  a  researcher  can  calculate  the  transmitted  and  reflected  amplitudes  of  the  electromagnetic 
fields  of  a  beam  incident  on  the  stack,  if  the  film  thicknesses  and  the  complex  indices  of  each  film 
and  those  of  the  incident  and  exit  media  are  known.  From  these  amplitudes,  one  is  able  to  find  the 
absolute  transmittance  and  reflectance  of  the  beam  for  either  mode  and  for  an  unpolarized  iwgm 
(Appendix  B). 

We  distinguish  two  cases  of  reflectance  and  transmittance  depending  on  whether  or  not 
interference  effects  occur  in  the  substrate. 

3.2  OPTICAL  MODEL  IN  THE  CASE  OF  PHASE  INCOHERENCY  THROUGHOUT 
THE STACK 


The  transfer  matrices  applicable  to  a  beam  traversing  a  stack  in  which  there  is  no  interference 
relate  beam  intensities  at  various  places  in  the  stack  and  contain  no  complex  numbers.  We  must, 
however,  continue  to  distinguish  between  the  modes  of  the  incident  beam  as  in  Sections  3.1.1  and 
3.1.2  above  (Appendix  C). 

3.2.1  Phase  Incoherent  Transfer  Matrix  at  an  Interface  of  Two  Media 

Suppose  l/  represents  the  beam  intensity  in  the  jth  medium  going  in  the  +z  direction  at 
the  interface  of  the  jth  and  mth  media,  and  Ij'  is  the  beam  intensity  in  the  -z  direction.  If 

corresponding  intensities  of  the  beam  in  the  mth  medium,  the  transfer  matrix 
among  these  intensities  is  (Appendix  C) 


-jm 


c 


where  the  intensity  transfer  matrix  Sj„  is 


-jm 


T. 

JD 


~R. 


ml 


L^Jm  ~ 


(18) 
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The  traosmittaace  (Tj„)  and  reflectance  (Rj„)  coefficients  come  from  the  Fresnel  ampli¬ 
tude  coefficients  (Eq.  (16)), 


T,.  =  KUr'J] 


(19) 


where  Re  indicates  the  real  part  of  a  complex  number  and 

^jin  ■" 


3.2  J  Phase  Incoherent  Transfer  Matrix  Through  a  Uniform  Medium 

The  phase  incoherent  transfer  matrix  ^  between  the  intensities,  and  Ig'at  the  entrance 
surface  of  a  uniform  layer  and  the  corresponding  intensities,  and  I  ,  at  the  layer’s  exit  surface 
is  (Appendix  C) 


and 


= ; 


c  = 


exp  (2Im  [6] )  0 

0  exp  (-2Im  [5]  )J 


(20) 


Im  refers  to  the  imaginary  part  of  a  complex  number. 

3.23  The  Transmittance  and  Reflectance  of  a  Beam  Incident  on  Thin  Films  Deposited 
on  a  Substrate  in  the  Case  of  Phase  Incoherency  Throughout 

The  transmittance  T  and  reflectance  R  of  a  beam  incident  upon  a  layered  system  of  films  on 
a  substrate  come  directly  from  the  intensity  transfer  matrices  and  ^  in  the  case  there  is  inter¬ 
ference  in  any  film  or  the  substrate  (Appendix  C).  As  in  the  phase  coherent  case,  the  transmittance 
and  reflectance  of  waves  in  the  TE  mode  differ  from  those  in  the  TM  mode. 


33  COMPUTER  PROGRAM  STACK 


The  model  of  the  transmittance  T  and  reflectance  R  of  a  light  beam  incident  on  a  stack  of 
uniform,  homogeneous,  layered  materials  consisting  of  films  on  either,  or  both,  sides  of  a  substrate 
is  used  in  the  computer  program  STACK  to  compute  T  and  R.  To  use  STACK,  one  needs  the 
incident  beam  angle  and  the  indices  and  the  thicknesses  of  each  film  of  the  layered  system  and 
the  substrate. 
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STACK  users  can  select  one  of  three  coherence  cases: 

1.  Interference  occurs  in  all  films  and  the  substrate. 

2.  Interference  occurs  in  all  films  but  not  in  the  substrate. 

3.  No  interference  takes  place  in  any  film  or  the  substrate. 

STACK  is  written  in  the  Interactive  Data  Language  (IDL)  and  is  a  replacement  for  program 
CALCRT  (AEDC  Technical  Memorandum  Report  of  May  1986),  which  calculates  T  and  R  for  a 
beam  incident  on  a  single  film  on  a  substrate  (See  Appendix  D).  The  model  of  a  single  film  on  a 
substrate  is  also  used  in  AEDC  program  TRNLIN  (Ref.  7)  which  finds  the  optical  indices  of  a  film 
from  normal  transmittance  measurements  of  an  incident  beam.  To  date  no  AEDC  programs  exist 
which  compute,  from  transmittance  measurements,  the  optical  indices  for  a  film  which  is  stacked 
with  other  films  on  a  substrate. 

4.0  EXAMPLES  OF  ANALYSES  OF  OPTICAL  DATA  OF  MIXTURES 

For  the  past  several  years,  AEDC  researchers  have  measured  the  infrared  optical  constants 
of  many  materials  condensed  on  ciyogenically  cooled  substrates.  Some  of  the  data  are  for  the  pure 
form  of  materials  at  different  temperatures  (and  densities),  and  some  are  measurements  of  the  op¬ 
tical  properties  of  mixtures  of  these  pure  materials  (Refs.  2-5,  8,  and  9).  We  have  used  these  data 
to  illustrate  the  Lorentz-Lorenz  transformation  of  the  optical  constants  of  a  single  material  and 
some  analysis  methods  for  the  optical  properties  of  mixtures  of  materials. 

4,1  LORENTZ-LORENZ  TRANSFORMATIONS  OF  n  OF  A  GIVEN  MATERIAL  AT 
DIFFERENT  DENSITIES 

The  Lorentz-Lorenz  constant  of  Eq.  (5)  allows  researchers  to  compute  the  optical  constants 
of  a  pure  material  at  one  density  from  their  values  at  another  density.  We  applied  the  Lorentz- 
Lorenz  transformation  to  the  optical  constants  measured  at  AEDC  for  three  pure  materials  (COj. 
HjO,  and  NH3)  which  had  formed  cryofilms  on  a  gennaniurn  (Ge)  substrate  at  two  temperatures, 
20K  and  80K.  Consequently,  we  could  calculate  the  optical  constants  of  each  material  at  80K  from 
their  values  at  20K  and  check  the  calculated  values  against  the  measured  80K  values.  We,  of 
course,  could  go  the  other  way  and  compute  20K  values  of  the  indices  from  the  measured  80K  val¬ 
ues  and  compare  with  those  actually  observed  at  20K. 
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4.1.1  The  Lorentz-Lorenz  TransfonnatioD 


The  transfonnation  of  the  complex  optical  index  n  of  a  material  with  a  (mass)  density  r  to 
the  optical  index  n'  that  the  material  has  when  its  density  is  r'  comes  from  the  fact  that  the 
Lorentz-Lorenz  constant  (Eq.  (5))  must  not  change  when  the  primed  quantities  replace  the 
unprimed  ones.  We  assume  that  the  local  field  parameter  b  does  not  differ  with  density.  Equating 
the  primed  and  unprimed  Lorentz-Lorenz  constants,  one  can  find  n’  from 


fi-^be: 

p 

J  \  p  / 

nM-P+b-HP- 
P  P 


(21) 


The  computer  program  LLTRNSFM,  written  in  IDL,  uses  Eq.  (21)  to  compute  the  n  and  k 
spectra  of  a  material  at  a  given  density  from  the  measured  n  and  k  spectra  at  a  different  density. 


4.1.2  Lorentz-Lorenz  Transformations  of  the  Optical  Constants  of  COj  Films 

Researchers  at  AEDC  have  observed  the  infrared  n  and  k  spectra  of  CO2,  H2O.  and  NHj  films 
formed  on  a  Ge  substrate  at  20K  and  also  at  80K.  They  used  program  TRNLIN  to  find  least-squares 
values  of  the  optical  constants  from  normal  transmittance  measurements  and  then  converted  the 
least-squares  k  spectrum  to  an  n  spectrum  via  the  subtractive  Kiamers-Kronig  transformation  in 
program  SKKK2N. 


n(v)  =  n(v„)  + 


k(v')-vk(v)  v'k(v') -v„k(v„)- 


A  2 
V  -V 


,2  2 

V  -v„ 

m 


dv'. 


(22) 


is  a  reference  wave  number,  and  P  denotes  the  Cauchy  principal  integral  value. 

4.1.2.1  Transformation  of  n  of  CO2  from  20K  to  80K 

The  authors  transformed  the  20K  Kramers -Kronig  n  values  of  CO2  and  its  least-squares  k 
values  to  those  n  and  k  spectra  that  the  Lorentz-Lorenz  transformation  predicts  at  80K.  The  density 
of  the  CO2  films  was  1080  kg/m^  at  20K  and  1670  kg/m^at  80K  (Table  1). 

The  transformed  (Lorentz-Lorenz)  n  spectrum  (for  b  =  2)  near  the  strong  'HTOj  at  2350 
cm'’  appears  in  Fig.  4  along  with  the  observed  (Kramers-Kronig)  values  of  n  observed  at  80K; 
the  analogous  transformed  and  observed  (least-squares)  k  spectra  are  in  Fig.  5.  Away  from  strong 
absorption  bands,  the  transformed  and  observed  n  values  are  in  good  agreement.  However,  the 
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peak  value  of  k  (Loientz)  is  shifted  downward  in  wave  number  from  the  measured  k  peak,  and  the 
extrema  of  n  (Lorentz)  are  displaced  to  lower  wave  numbers  from  those  observed  for  n  at  80K.  In 
addition,  the  peak  of  the  Lorentz-Lorenz  value  of  k  is  greater  than  its  measured  value,  and  the  de¬ 
viations  in  the  Lorentz-Lorenz  n  spectrum  are  more  than  the  deviations  observed  for  n  at  80K. 

To  check  the  effect  of  measurement  errors  in  the  density  values  of  the  COj  films,  we  recal¬ 
culated  the  Lorentz-Lorenz  spectra  with  new  density  values  and  noted  that  as  the  ratio  of  the 
densities,  p  (80K)/  p  (20K),  decreases  towards  unity,  the  spectral  shifts  lessened  but  the  agreement 
between  the  n  spectra  away  from  the  strong  band  worsened.  We  conclude  from  this,  and  other 
examples,  that  the  better  the  n  values  (away  from  strong  absorptions)  agree,  the  more  confidence 
we  have  that  the  density  values  are  correct. 

We  also  experimented  with  the  value  of  the  local  field  parameter  b  and  found,  to  our  surprise, 
that  if  just  the  value  of  b  is  changed  from  2  to  10™  (essentially  infinity),  the  n  and  k  spectra  resulting 
from  the  transformation  agree  surprisingly  well  with  the  corresponding  observed  spectra  (Figs.  6 
and  7).  This  means,  by  Eq.  (2),  that  parameter  a  in  Eq.  (3)  is  zero  in  this  example.  Thus,  in  contrast 
to  the  normal  assumption  for  non-polar  lih  dielectrics  (where  a  -  S/Cq  and  b  -  2),  the  present  result 
indicates  that  the  (average)  contribution  of  the  polarization  vector  P  to  the  local  electric  field  Ei„c 
within  the  dielectric  is  negligible  for  CO2  films. 

4.12.2  Transfonnatioa  of  n  of  COj  firom  80K  to  20K 

The  n  and  k  spectra  that  come  from  a  Lorentz-Lorenz  transformation  of  80K  n  and  k  data  to 
their  values  at  20K  are  pictured  in  Figs.  8  and  9.  In  contrast  to  the  20K  to  80K  transformation,  the 
Lorentz-Lorenz  spectra  have  shifts  upward  in  wave  number,  and  the  peak  values  of  k  and  n  are 
smaller  than  those  experimentally  observed.  As  in  the  20K  to  80K  transformation,  the  spectral 
shifts  decrease  as  the  ratio  of  the  densities  decreases  toward  one,  but  the  n  values  away  from  the 
strong  absorption  band  diverge.  Again,  when  the  value  of  b  is  10™,  the  corresponding  transformed 
and  observed  spectra  become  more  alike  (Figs.  10  and  11). 

4.13  Lorentz-Lorenz  Transformations  of  the  Optical  Constants  of  H^O  Films 

We  applied  the  above  analysis  techniques  for  COj  to  20K  and  80K  n  and  k  spectra  of  H2O 
films.  The  density  of  20K  film  (670  kg/m^)  is  agmn  greater  than  the  80K  film  density  (920  kg/m^) 
(Table  1). 
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4.13.1  Transforniation  of  n  of  HjO  fh>in  20K  to  80K 

The  Lxtrentz-Lorenz  transformed  n  and  k.  spectra  (from  20K  to  80K)  are  compared  to  the 
corresponding  observed  values  in  Figs.  12  and  13.  Though  the  two  k  spectra  are  close  in  the  strong 
absorption  band  centered  near  3300  cm'‘,  they  differ  greatly  in  the  1200  to  2800  cm'’  and  750  cm'* 
regions,  and  the  transformed  n  spectrum  is  everywhere  greater  than  the  observed  one.  By  adjusting 
the  values  of  either  the  densities  or  of  the  parameter  b,  one  can  get  closer  agreement  in  certain 
regions  of  the  spectra  but  at  the  expense  of  worse  agreement  in  other  regions.  The  cause  for  this 
may  be  that  the  physical  structure  of  the  HjO  films  is  different  at  the  two  temperatures  much  like 
the  case  of  NHj  films.  (See  Section  4.1 .4  below.)  It  is  also  possible  that  the  analysis  methods  need 
revision  or  the  original  measurement  and  data  reduction  procedures  are  in  error. 

4.13.2  Transformation  of  n  of  HjO  from  80K  to  20K 

Figures  14  and  15  compare  the  transformed  n  and  k  spectra  to  the  observed  ones  when  going 
from  80K  to  20K.  Again,  the  agreement  between  the  k  spectra  is  best  near  3300  cm'’,  but  poor  else¬ 
where,  and  the  n  spectra  are  not  close  anywhere.  No  value  of  the  parameter  b  allows  better  agree¬ 
ment  everywhere  in  the  spectra;  the  same  is  true  for  the  densities. 

4.1.4  Lorentz-Lorenz  Transformations  of  the  Optical  Constants  of  NH3  Films 

The  Lorentz-Lorenz  n  and  k  spectra  and  the  corresponding  observed  spectra  for  NH3  films 
are  in  Figs.  16  and  17  for  the  20K  to  80K  transformation  and  in  Figs.  18  and  19  for  the  80K  to  20K 
transformation.  The  film  density  at  20K  (670  kg/m^)  is,  as  for  the  two  film  materials  above,  less 
than  that  at  80K  (920  kg/m^). 

Because  the  absorption  bands  measured  for  the  80K  films  are  much  stronger  and  sharper 
than  those  measured  at  20iC,  the  transformed  k  spectra  bear  little  resemblance  to  the  observed  k 
spectra.  However,  except  near  the  bands  at  1 100  cm  ’  and  3350  cm'’,  the  n  spectra  agree  well. 
Varying  b  or  the  density  values  does  not  change  the  computed  k  spectra  enough  to  mimic  the 
observed  k  spectra. 

The  major  reason  for  this  nonagreement  between  the  20K  and  80K  NHj  spectra  is  that  the 
structure  of  the  80K  film  was  crystalline  while  that  of  the  20K  film  was  amorphous  or  metastable 
(Ref.  10).  This  accounts  for  the  stronger  and  narrower  features  in  the  experimental  80K  spectra 
compared  to  those  in  the  20K  spectra;  the  Lorentz-Lorenz  transformation  does  not  change 
significantly  the  sharpness  of  the  peaks. 
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4.1.5  Discussion  of  Lorentz-Lorenz  Transformation  - 

The  results  of  the  Lorentz-Lorenz  transformation  for  COj  films  show  that  this  analysis 
method  may  yield  n  and  k  spectra  similar  to  those  actually  measured  at  a  different  density  provided 
the  film  material  is  predominantly  made  of  molecules  which  do  not  form  hydrogen  bonds,  and  the 
proper  value  of  the  local  field  parameter  b  is  used.  For  films  composed  of  hydrogen-bonding 
molecules  like  HjO  and  NHj,  the  physics  is  more  complicated  than  the  simple  changes  of  the  optical 
constants  that  density  differences  alone  can  create.  To  create  an  adequate  model,  one  probably  must 
look  in  (quantum  mechanical)  detail  at  the  microstructure  of  films  of  hydrogen-bonding  molecules 
and  the  changes  that  occur  because  of  different  condensation  temperatures,  molecular  migration 
rates,  annealing  effects,  and  so  forth. 

4.2  LORENTZ-LORENZ  ANALYSIS  METHODS  OF  THE  OPTICAL 
PROPERTIES  OF  MIXTURES  OF  MATERIALS 

One  of  the  goals  of  the  AEDC  research  program  on  contaminants  of  spacecraft  is  to  predict 
the  optical  properties  of  a  film  composed  of  a  mixture  of  materials  from  the  known  optical  indices 
of  each  component  material.  In  this  section,  we  calculate  optical  constants  of  several  cryogenic 
films  which  are  mixtures  of  some  simple  molecules,  with  the  help  of  the  Lorentz-Lorenz  result  of 
Eqs.  (8)  and  (9),  and  compare  these  indices  with  the  n  and  k  spectra  measured  for  the  mixture. 

We  also  compared  the  computed  transmittance  spectra  of  selected  films  of  mixtures  with 
their  computed  transmittance  spectra  calculated  with  measured  optical  constants  of  the  mixture  and 
some  spectra  computed  as  if  the  components  of  the  mixture  formed  separate  layers  within  the  larger 
film.  To  compute  the  last  spectra,  we  must  decide  what  thickness  value  each  layer  might  have  arid 
in  what  order  are  the  layers  of  the  film,  none  of  which  comes  from  experimental  data.  We  make  the 
assumption  that  each  density,  p  ^  that  component  i  has  in  its  layer  of  the  film  is  related  to  its  density 
p  pi  in  its  pure  state  by  the  same  factor  as  any  other  component;  that  is,  pi  =  pppj ,  where  P  is  con¬ 
stant  for  all  components  L  When  we  combine  this  assumption  with  the  fact  that  tte  sum  of  the  thick¬ 
ness  values,  d,,  of  each  layer  must  be  the  (measured)  value  of  the  total  film  thickness  d,  we  have 

d/d  =  (M/M)  (p/(pppi)), 


where 


P  =  (S(M/M)  (p/ppi)  ).  (23) 

Here.  Mj/M  is  the  mass  ratio  of  the  ith  component  and  p  is  average  density,  respectively,  of  the 
film.  (If  the  mole  fractions,  Ni/N,  are  known  instead  of  mass  ratios,  we  can  use  the  conversion. 


19 


AEDC-TR-94-3 


=  (N/N)  (m/m),  where  in/m  is  the  molecular  mass  ratio  of  the  ith  component  to  the  average 
molecular  mass  of  the  components.) 

We  also  assumed  that  the  film  component  whose  pure  state  equilibrium  vapor  pressure  was 
the  lowest,  at  a  given  substrate  temperature,  would  form  the  first  layer  on  the  substrate;  the 
component  with  the  next  lowest  vapor  pressure  would  form  the  next  layer  on  top  of  the  first  layer, 
and  so  forth.  Thus,  by  the  above  criterion,  for  a  Nj/HjO/COj/CO  film  structured  in  layers,  the  order 
of  the  layers  (starting  with  the  side  on  which  the  beam  is  incident)  is  CO,  Nj,  CO2,  and  HjO. 

4.2.1  Optical  Properties  of  20K  50-Percent  COi/SO-Percent  CO  Films 

Several  years  ago,  researchers  at  AEDC  measured  the  optical  properties  of  films  made  of 
small  molecules  condensed  onto  a  Ge  substrate  at  20K  or  80K  (Ref.  5).  We  used  these  data  in 
program  NMIX  to  compute  the  optical  constants  of  20K  films  composed  of  SO-peicent  (mole 
fraction)  COj  molecules  and  50-percent  CO  molecules  and  to  compare  to  the  measured  n  and  k 
spectra;  all  b  values  were  set  to  2.  The  measured  and  computed  n  and  k  spectra  appear,  respectively, 
in  Figs.  20  and  21.  Except  near  the  large  CO2  absorption  band  at  2350  cm'’,  the  calculated  and 
measured  n  spectra  are  close,  though  usually  the  calculated  n  values  are  greater  than  the  correspond¬ 
ing  measured  values.  Both  the  n  and  k  spectra  are  shifted  slightly  from  one  another  near  2350  cm''. 

The  measured  k  spectrum  is  smaller  than  the  corresponding  calculated  one  at  the  large  COj 
band  but  greater  at  the  2150  cm*'  CO  band.  This  may  mean  that  the  concentration  of  CO2  molecules 
is  less  than  50  percent.  (However,  when  the  mole  fraction  of  COj  is  set  to  40  percent,  so  that  the 
new  calculated  k  spectrum  nearly  matches  the  measured  spectrum  near  the  CO2  band,  at  the  CO 
band  near  2150  cm''  the  measured  k  values  are  still  much  greater  than  the  calculated  ones.)  Like¬ 
wise,  making  small  changes  in  the  values  of  b  improves  the  agreement  of  the  calcidated  n  and  k 
spectra  with  the  measured  spectra,  but,  in  contrast  to  the  case  of  pure  COj  films  discussed  above, 
when  b  is  set  to  10™,  the  calculated  spectra  are  highly  distorted  from  the  measured  ones. 

Figure  22  is  a  plot  of  three  computed  normal  transmittance  infirared  spectra  for  a  I3.16-)im 
CO2/CO  on  a  20K  Ge  substrate.  The  spectrum  labeled  ‘T  (obs)”  is  calculated  from  the  measured 
optical  constants  of  the  mixture  (Ref.  5);  ‘T  (mix)”  is  computed  from  the  calculated  optical 
constants  (from  program  NMIX)  of  the  mixture;  and  ‘T  (layers)”  is  a  STACK  calculation  of  the 
transmittance  as  if  a  uniform  6.08-pm  CO  layer  were  next  to  a  7.08-)im  layer  of  COj  which  is  on 
the  substrate.  T  (mix)  and  T  (layers)  are  nearly  identical,  which  may  seem  surprising  because  in 
the  latter  case,  there  are  actually  two  films  stacked  on  the  substrate,  and  multiple  reflections  can 
occur  at  their  mutual  interface.  Nevertheless,  the  interference  fringes  of  these  two  spectra  mimic 
one  another.  This  is  probably  because  the  refractive  indices  of  each  of  the  layers  are  close  to  each 
other  in  value  and  that  of  Ge  is  several  times  greater.  The  period  of  the  T  (obs)  spectrum  is  notice- 
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ably  greater  than  those  of  T  (mix)  and  T  (layers)  and  indicates  the  refractive  index  of  the  observed 
film  is  less  than  that  calculated  in  NMIX.  These  spectra  are  very  similar  to  the  measured  spectrum 
found  in  Fig.  1 1  of  Ref.  S. 

4  J  J  Optical  Properties  of  20K  64-Peiiceiit  Nj/ZS-Percent  CO/13-Percent  COj  Films 

When  Nj,  which  has  no  near  infrared  absoiption  bands,  is  added  to  the  mixture  of  CO  and 
CO2,  the  expected  n,  k,  and  T  spectra  are  similar  to  those  of  the  CO2/CO  films.  Indeed  this  is  die 
case,  as  Figs.  23,  24,  and  25  show  for  n,  k,  and  T,  respectively. 

Away  from  the  COj  band  at  2350  cm'\  the  calculated  n  values  roughly  agree  with  the 
measured  values,  though  they  are  greater  over  most  of  the  spectrum.  The  extrema  of  the  measured 
n  spectrum  are  much  more  pronounced  near  the  COj  band  and  are  shifted  from  the  calculated  val¬ 
ues  as  well.  The  measured  k  spectrum  has  peak  values  greater  than  the  calculated  ones  at  the  2150 
cm*’  CO  band  and  at  2350  cm*'.  The  differences  in  the  calculated  n  and  k  spectra  fiom  the  mea¬ 
sured  spectra  indicate  that  the  concentrations  of  both  CO  and  CO2  molecules  may  have  been  higher 
than  23  percent  and  13  percent,  respectively.  (The  mole  fractions  were  found  from  chemical  anal¬ 
ysis  of  the  gaseous  mixture  used  to  make  the  film  deposit.  The  actual  composition  of  the  film  on 
the  substrate  cannot  be  measured  directly.) 

In  Hg.  25,  the  T  (layers)  spectrum  is  calculated  in  STACK  as  if  the  6.602-pm  film  had  a  uni¬ 
form  CO  layer  1.521  pm  thick  at  its  outer  surface  (closest  to  the  source  of  the  infrared  beam),  a 
1.001-)im  layer  of  CO2  next  to  the  substrate,  and  a  4.080-pm  Nj  layer  sandwiched  between  them. 
Again,  it  is  noteworthy  that  though  interference  effects  can  occur  within  each  layer,  the  interference 
fringe  pattern  has  a  structure  and  a  period  that  is  very  much  like  the  one  in  the  T  (obs)  spectrum  in 
which  the  observed  film  optical  constants  were  used  and  the  deposit  was  treated  as  a  single  film 
The  single  film  calculation  using  the  optical  constants  from  NMDC  produced  a  spectrum ‘T  (mix)" 
which  has  a  fringe  period  less  than  those  of  the  other  two.  This  indicates  the  values  of  n  from  NMIX 
are  greater  than  the  observed  values  and  may  be  the  result  of  possible  errors  in  the  concentration 
values  of  the  components  mentioned  above.  The  measured  T  spectrum  of  the  film  (Fig.  13  of  Ref. 
5)  has  features  that  are  most  like  those  of  T  (obs). 

4,2.3  Optical  Properties  of  20K  S03>Percent  N2/22.5-Percent  H20/17,2-Percent  CO2/IO.O- 
Percent  CO  Films 

Our  last  example  of  a  20K  film  is  a  mixture  of  the  components  of  the  last  film  material,  N2, 
CO2,  and  CO,  with  HjO  molecules  that  can  form  hydrogen  bonds  with  other  molecules.  (This 
mixture  is  named  “Simulated  Plume  Mixture  1”  in  Ref.  5.)  The  n,  k,  and  T  spectra  of  these  films 
appear,  respectively,  in  Figs.  26, 27,  and  28. 
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In  most  Spectral  regions  away  from  the  COj  band  at  2350  cm  ',  the  calculated  n  and 
measured  n  spectra  are  in  good  agreement.  Near  2350  cm  ',  the  peaks  of  the  calculated  k  spectrum 
are  much  less  than  those  measured,  and  the  calculated  peak  is  shifted  upwards  in  wave  number  from 
2346  to  2360  cm"'.  However,  the  calculated  peak  k  values  at  the  water  band  near  3300  cm"'  are 
much  greater  than  the  measured  peaks.  There  are  at  least  two  causes  for  these  discrepancies 
between  spectra.  First,  the  concentration  values  of  the  absorbing  molecules  may  be  in  error.  (It 
appears  that  the  concentration  of  CX3j  may  be  higher  than  17.2  percent,  and  that  of  HjO  may  be 
lower  than  22.5  percent.)  Second,  and  probably  more  importantly,  the  degree  of  hydrogen  bonding 
of  the  water  molecules  is  highly  dependent  on  its  concentration  (Ref.  S).  As  the  concentration  of 
HjO  molecules  grows,  the  3300  cm''  0-H  stretch  band  rapidly  becomes  more  intense  and  wider  in 
a  nonlinear  way.  Evidently,  for  this  film  material  the  degree  of  hydrogen  bonding  is  much  less  than 
that  observed  in  pure  HjO  films  (Fig.  13). 

The  differences  between  the  calculated  and  observed  n  and  k  spectra  are  evident  in  the 
calculated  transmittance  spectra  of  a  20K  Nj/HjO/COjAX)  film  6.72-iim  thick  that  appear  in  Fig. 
28.  The  spectrum  calculated  from  observed  film  optical  constants,  T  (obs),  matches  the  measured 
spectrum  very  well  (Fig.  21  of  Ref.  5).  That  calculated  from  the  optical  constants  computed  in 
NMIX,  T  (mix),  is  similar  to  T  (layers).  (In  the  calculation  of  the  T  (layers)  spectrum,  the  order  of 
the  layers  is:  CO  (0.70  pm),  Nj  (3.40  pm),  COj  (1.40  pm),  HjO  (1.21  pm),  and  Ge  substrate.)  The 
T  (mix)  and  T  (layers)  spectra  have  stronger  and  broader  HjO  bands  than  observed  because  their 
calculations  are  based  on  the  optical  constants  of  pure  HjO  spectra  where  the  degree  of  hydrogen 
bonding  is  at  a  maximum. 

4,2,4  Optical  Properties  of  80K  Films  of  a  Mixture  of  Uralane^ 

andRTV®-566 

Our  last  example  is  a  mixture  of  47  percent  (by  weight)  Uralane  5753-A/B  (LV)  urethane 
casting  compound  and  53  percent  (by  weight)  RTV-566  silicon  sealant  (“Mixture  2”  in  Ref.  4),  The 
measured  and  calculated  n,  k,  and  T  spectra  are  in  Figs.  29-31. 

The  calculated  n  spectrum  is  larger  everywhere  than  the  measured  n  spectrum;  the  calmiat^ 
k  spectrum  is  greater  than  measured  k  spectrum  except  at  3300  cm  '  where  they  are  essentially 
identical.  The  T  (obs)spectium,  computed  from  the  observed  optical  constants,  is  in  good  agree¬ 
ment  with  the  measured  T  spectrum,  T  (meas),  though  near  2000  cm"'  T  (obs)  is  greater  than 
T(meas).  The  T  (mix)  spectrum  is  less  than  the  measured  T  spectrum,  but  they  agree  well  at  3300 
cm''.  The  agreement  of  the  spectra  can  be  improved  by  redoing  the  calculations  with  a  slight 
change  in  the  concentrations  of  the  components.  The  T  (layers)  spectrum  is  for  a  uniform  0.894- 
|ira  layer  of  RTV-566  on  a  0.692-pim  layer  that  is  next  to  the  80K  Ge  substrate.  Of  the  three 
calculated  spectra,  it  has  the  worse  agreement  with  the  measured  T  spectrum. 
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4  J.5  Discussion  of  the  Lorenti-Lorenz  Methods  for  Computing  the  Optical 
Properties  of  Mixtures 

The  Lorentz-Lorenz  methods  to  predict  the  optical  constants  of  a  film  composed  of  a 
mixture,  whose  optical  constants  in  the  pure  state  are  known,  can  provide  results  which  are  accurate 
to  within  10  percent,  or  so,  away  from  strong  absorption  bands.  One  difficulty  in  getting  more 
accurate  predictions  is  the  experimental  error  in  the  values  of  the  concentrations  of  each  film 
omponent.  Other  difficulties  include  the  theoretical  and  experimental  uncertainties  of  the  effects 
of  intermolecular  attractions.  All  difficulties  become  more  acute  as  the  number  of  constituents  in 
the  mixture  increases. 

The  problem  of  intermolecular  attractions  is  compounded  by  the  phenomenon  of  hydrogen 
bonding  which  affects  many  common  molecules  like  HjO  and  NH3.  Because  spectral  features 
change  drastically  as  the  degree  of  hydrogen  bonding  varies  (Ref.  11),  the  optical  constants  of  the 
pure  state  of  hydrogen-bonding  molecules  do  not  have  the  spectral  features  needed  for  materials 
where  the  degree  of  hydrogen  bonding  is  less.  To  predict  the  optical  constants  in  these  cases,  one 
might  use  the  optical  constants  measured  for  hydrogen-bonded  molecules  diluted  in  an  inert  matrix 
such  as  solid  Nj  or  Ar  For  instance,  a  researcher  might  be  able  to  use  the  r^tical  constants  of  a  film 
of  the  appropriate  concentration  of  H2O  molecules  in  N2  to  predict  more  accurately  the  spectral 
features  due  to  HjO  in  the  20K  Nj/HjO/COj/CO  film  (Figs.  26-28). 
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Figure  1.  Amplitudes  and  intensities  at  an  interface. 
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Figure  2.  Amplitudes  and  intensities  in  a  film  medium. 
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Figure  3.  Stack  of  films  on  a  substrate. 
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igure  4.  Lorentz-Lorenz  transformation  of  n  for  CO2  films  (20K  to  80K). 
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Figure  5.  Lorentz-Lorenz  transformation  of  k  for  CO2  films  (20K  to  80K). 
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Figure  6.  Lorentz-Lorenz  transformalion  of  n  for  COj  films 
(20K  to  80K,  b  -  1.0E20). 
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Figure  7.  Lorentz-Lorenz  transformation  of  k  for  CO2  films 
(20K  to  80K,  b  -  L0E20). 
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Figure  8.  Lorentz-Lorenz  transformation  of  n  for  CO2  films 
(80K  to  20K,  b  -  1.0E20). 
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Figure  9.  Lorentz-Lorenz  transformation  of  k  for  CO2  films 
(80K  to  20K,  b  -  1.0E20). 
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Figure  10.  Lorentz-Loreoz  transformation  of  n  for  CO2  films 
(80K  to  20K,  b  =  1.0E20). 
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Figure  11.  Lorentz-Lorenz  transformation  of  k  for  CO2  films 
(SOK  to  20K,  b  -  1.0E20). 
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Figure  12.  Lorentz-Lorenz  transformation  of  n  for  H2O  films 
(20K  to  80K). 


(20K  to  80K). 
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Figure  14.  Lorentz-Lorenz  transformation  of  n  for  H2O  films 
(80K  to  20K). 


Figure  15.  Lorentz-Lorenz  transformation  of  k  for  H2O  films 
(80K  to  20K). 
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Figure  16.  Comparison  of  n  (meas)  and  n  (Lorenz)  of  80K  NH3  films 
(from  20K  data). 


Figure  17.  Comparison  of  k  (meas)  and  k  (Lorenz)  of  80K  NH3  films 
(from  20K  data). 
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Figure  18.  Comparison  of  n  (meas)  and  n  (Lorenz)  of  20K  NH3  films 
(from  80K  data). 
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Figure  19.  Comparison  of  k  (meas)  and  k  (Lorenz)  of  20K  NH3  films 
(from  80K  data). 
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Figure  23.  n  (calc)  and  n  (meas)  of  20K  N2/CO/CO2  film. 
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Figure  24.  k  (calc)  and  k  (meas)  of  20K  N2/CO/CO2  film. 
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Figure  25.  Calculated  transmittance  of  6.602-pm  N2/CO/CO2  film. 
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Figure  26.  n  (calc)  and  n  (meas)  of  20K  N2/H2O/CO2/CO  film. 


0  1000  2000  3000 

Wavenumber,  cm'' 

Figure  27.  k  (calc)  and  k  (meas)  or20K  N2/H2O/CO2/CO  film. 
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Figure  28.  Calculated  transmittance  of  6.72 -Itm  N2/H2O/CO2/CO  film. 


Figure  29.  n  (calc)  and  n  (meas)  of  80K  mixture  2  (Uralane®  and  RTV®)  films. 
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Figure  31.  Measured  and  calculated  transmittance  of  80K  mixture  2 
(Uralane®  and  RTV®)  films. 
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Ikble  1.  Densities  of  Selected  FUm  Materials 


Film 

Ibrnperatuie 

Density  (kg/m^) 

CO2 

20K 

1080 

CO2 

80K 

1670 

H2O 

20K 

670 

H2O 

80K 

920 

NH3 

20K 

760 

NH3 

80K 

870 

N2 

20K 

830 

CO 

20K 

800 

Uralane® 

80K 

863 

RTV®.566 

80K 

754 

coyco 

20K 

950 

N2/CO/CO2 

20K 

941 

20K 

770 

Mixture  2 

80K 

938 
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APPENDIX  A 

The  Lorentz-Lorenz  Rebtion 

This  appendix  contains  an  outline  of  the  derivation  of  the  classical  LoFentz-Lorenz  relation 
which  relates  the  concentration  of  molecules  of  a  Uh  dielectric  to  its  complex  refractive  index.  The 
discussion  is  similar  to  those  given  by  Guenther  (Ref.  12)  and  Bom  and  Wolf  (Ref.  13). 

Suppose,  for  the  moment,  that  each  molecule  in  the  dielectric  has  a  single  charge,  q,  of  mass, 
|i,  which  acts  as  a  damped  harmonic  oscillator  to  give  the  optical  properties  of  the  dielectric.  If  x 
is  the  displacement  of  the  oscillator  from  equilibrium  and  E|g,.  is  the  component  of  the  local  electric 
field  experienced  by  the  oscillator,  Newton’s  Second  Law  gives  the  following  equation  of  motion: 

X  +  T*  =  tojx  =  ^,  (A-1) 

where  y  and  tOg  are  the  damping  and  simple  harmonic  frequency  constants,  respectively. 


To  make  the  connection  between  macroscopic  and  microscopic  quantities,  one  can  use  the 
polarization  density  P  s  Np,  where  N  is  the  number  density  of  the  oscillators  and  p  =  qx  is  the 
electric  dipole  moment  per  oscillator.  The  local  electric  field  is  related  to  P  and  the  “cavity”  (or 
measured)  value  of  the  electric  field  E  by  =  E  +  aP  (Eq.  (3)).  Furthermore,  for  linear  dielectrics, 
E  and  P  are  linearly  proportional:  P  =  XE.  These  relations  allow  the  equation  of  motion  to  be  re¬ 
written  in  terms  of  P  and  its  time  derivatives.  If  we  make  the  assumption  that  the  incident  electric 
field  E  has  a  frequency  w  and  a  time  dependence  of  exp[-imt],  the  constant  c  can  be  found  fiom  the 
equation  of  motion  to  be 


X  = 


Ni 


tOo-to 


rNaq 


-It 


-lyto 


(A-2) 


The  relationsiop  of  X  to  die  complex  optical  index  n  of  a  dielectric  comes  from  the  “consti- 
tuitive”  relation  between  the  displacement  field  D  and  the  electric  field  E,  which  is  D  -  EgE  +  P 
=  eI,  where  e  is  the  (complex)  electric  permittivity  of  free  space.  This^ means  that  E  =  Eg  x  ■ 
Because  the  refractive  index  n  for  a  nonmagnetic  dielectric  is  n  ~  (e/Eg)2 ,  the  desired  relation  be¬ 
tween  x^dnis 


X  —  £  ■ 


(A-3) 
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When  one  equates  Eqs.  (A-2)  and  (A-3)  and  does  some  algebraic  manipulation,  the  result  is 


n^-1 
n*  +  b 


Nq^ 

EoH 


2  iNaq^ 
(Oo  -  (O'* — —  - 

-li 


(A-4) 


where  constant  b  is  given  by  Eq.  (4). 

Equation  (A-4)  is  one  form  of  the  Lorentz-Lorenz  relation,  but  the  relation  is  often  given  in 
terms  of  the  mean  polarizability  which  relates  the  (average)  electric  dipole  moment  per  molecule, 
p,  with  the  local  electric  field  Ei^c  in  a  lih  dielectric;  p  =  aE,^ ,  If  the  original  equation  of  motion, 
Eq.  (A-1),  is  recast  in  terms  ofa  and  the  time  derivatives  of  P,  a  turns  out  to  be 


a  = 


9! 

<o2  -  to^  -  iyoa 


(A-5) 


The  Lorentz-Lorenz  formula  of  Eq.  (1)  follows  when  Eq.  (A-5)  is  multiplied  by  Na  and  compared 
to  Eq.  (A-4), 


If  each  molecule  has  more  than  one  oscillator,  the  Lorentz-Lorentz  relation  of  Eq.  (1)  is  still 
valid  because  of  the  linear  relationship  of  the  electric  dipole  moment  p  and  the  local  electric  field 
El,,,  (See  Ref.  12). 


If  the  dielectric  is  a  homogeneous  mixture  of  non-inmracting  molecule  species,  each  with  a 
number  density  Nj  within  the  mixture,  the  mean  polarizability  of  the  mixture  a  is  the  mathematical 
weighted  mean  of  the  individual  polarizabilities  for  each  species  ,  a  =  (ZNjOj)  /N,  where  N  is 
the  total  number  density  of  molecules  in  the  mixture,  N  =  £Nj.  A  value  of  a^.  can  be  calculated 
from  Eq.  (1)  using  the  value  of  the  complex  index  nj  and  number  density  Npj  measured  for  the  jth 
molecular  species  j'n  its  pure  state.  When  these  relations  are  combined  with  the  corresponding 
equation  for  a  of  the  mixture,  a  is 


a  = 


nf  +hyj' 


(A-6) 


where  the  local  field  parameters  are  aj  and  bj  =  (£(,0^)  -  1 ,  for  the  jth  component  in  its  pure 
state,  and  a  and  b  =  (e^a)  '  -  1  for  tiie  mixture.  When  one  solves  Eq.  (A-6)  for  n,  Eqs.  (8)  and 
(9)  result. 
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APPENDIX  B 

Transmittance  and  Reflectance  of  a  Stack  of  Multilayers  In  the  Case  of  Phase 

Coherence  Within  the  Films 


Potter  (Ref.  6)  discusses  the  transmittance  and  reflectance  of  a  stack  of  uniform  multilayers 
on  a  substrate  (Fig.  3)  using  the  (arbitrary)  choice  of  exp[iwt]  for  the  time  dependence  of  an  incident 
light  beam.  This  forces  one  to  define  the  complex  refoactive  index  to  be  n  =  n  -  ik  and  makes  a 
difference  in  some  of  the  other  defined  quantities  and  derived  equations.  This  appendix  gives  the 
major  changes  in  the  article  by  Potter  that  result  from  a  field  time  dependence  choice  of  exp[-iwt] 
which,  among  other  things,  makes  n  =  n  +  ik.  The  stack  (in  Fig.  3)  has  Nj  uniform  films  ahead  of 
the  substrate  (labeled  s)  and  N,,  films  on  the  back  side  of  the  substrate. 

The  wave  equations  of  the  fields  that  come  from  Maxwell’s  Equations  do  not  change  with 
the  choice  of  time  dependence,  but  the  signs  of  the  exponents  for  the  field  quandties  in  the  TE  and 
TM  modes  are  reversed  as  Eqs.  (1()'14)  show  when  compared  to  the  corresponding  equations  in 
Ref.  6.  There  are  sign  reversals  as  well  in  Eqs.  (5a)  and  (5b)  of  Ref.  6  for  U’  and  V’  and  in  the  off- 
diagonal  terms  of  the  matrix  in  Eq.  (10)  of  Ref.  6.  The  transfer  matrix  (Eq.  (15))  is  unchanged 
from  Eq.  (13a)  of  Ref.  6,  but  the  exponents  of  L  (Eq.  (17))  have  their  signs  reversed  from  those  in 
Eq.  (14)  of  Potter  (Ref.  6). 

In  the  case  of  interference  within  the  films  and  substrate,  the  algorithm  in  program  STACK 
incorporates  Eqs.  (26),  (27),  and  supporting  equadons  of  Ref.  6;  if  interference  occurs  in  the  films 
but  not  the  substrate,  STACK  uses  Eqs.  (29)  and  (31)  of  Ref.  6  (and  their  supporting  equations). 
There  are  several  corrections,  however.  The  exponentials  in  Potter’s  Eq.  (30)  should  be  deleted, 
and,  most  importandy,  which  appears  in  Ref.  6  just  before  Eq.  (31),  has  the  definition, 
^/i  ■  represents  the  Fresnel  transmittance  coefficient  at  the  jth  interface  for 

a  beam  passing  from  the  sutetrate  to  the  front  surface  of  the  stack. 

Another  important  note  is  that  the  quantities,  ,  appearing  in  the  definition  of  the  transmis¬ 
sion  and  reflection  coefficients  (Eq.  (16))  depend  on  the  polarization  mode  of  the  beam.  Because 
only  ratios  of  appear  in  the  transmittance  and  reflectance  equations,  multiplicative  constants 
common  to  all  of  them  can  be  neglected.  In  the  case  of  nonmagnetic  materials,  is  proportional 
to  \/yj  in  the  perpendicular  (TE)  mode  and  proportional  to  -DjV  Yj  in  the  parallel  (TM)  mode.  Yj 
can  be  found  from 

(noS^nW^  (B-1) 

where  iiq  is  the  (real)  refiaction  index,  and  <|)o  is  the  beam  angle  in  the  incident  medium,  and  is 
the  complex  refraction  index  of  the  jth  medium. 
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APPENDIX  C 

Transmittance  and  Reflectance  of  A  Stack  of  Multilayers  in  the  Case  of  No 
Phase  Coherence  Within  the  Films  or  Substrate 


This  appendix  outlines  an  extension  of  Potter’s  analysis  to  the  case  of  a  light  beam  incident 
on  films  supported  by  a  substrate  in  which  no  interference  occuned  anywhere  in  the  stack.  In  this 
case,  the  transfer  matrices  relate  beam  intensities  at  various  points  within  the  stack  instead  of  beam 
amplitudes,  as  in  the  case  of  film  phase  coherence  (Figs.  1-3).  As  above,  the  stack  has  uniform 
films  ahead  of  the  substrate  (labeled  s)  and  N,,  films  on  the  back  side  of  the  substrate. 


For  either  polarization  (TEM  or  TMM),  the  reflectance  and  transmittance  of  a  beam  incident 
on  the  stack  are,  respectively. 


1  +  r.oRX  ’ 


(C-1) 


and 


To.t^e; 

1  +  r^rX 


(C-2) 


Rq,  and  Tg,  are  the  reflectance  and  transmittance,  respectively,  of  a  beam  going  fiom  the  non- 
absorbing  incident  medium  (semi-infinite  medium  0)  into  the  substrate,  R^  and  are  the 
corresponding  quantities  for  a  beam  originating  in  the  substrate  and  traveling  to  the  fiont  surface, 
and  R^  and  refer  to  analogous  quantities  of  a  beam  originating  in  the  substrate  and  entering  the 
exit  medium  (semi-infinite  medium  f).  In  the  jth  medium,  the  exponential  decay  term  Ej  is 

Ej=exp  [-2Im(5j)],  (C-3) 

where 


5j  =  27tvYj  d., 

Yj  isdefinedinEq.(B-l),  and  dj  is  the  thickness  of  the  jth  medium.  Im  refers  to  the  imaginary  part 
of  a  complex  number. 

Dq,  and  the  various  transmittances  and  reflectances  in  Eqs.  (C-1)  and  (C-2)  come  from  the 
intensity  transfer  matrices  for  the  stack.  These  matrices  are  themselves  products  of  the  intensity 
matrices  at  an  interface  and  through  a  uniform  medium  defined  in  Eqs.  (18-20).  The  total  transfer 
matrix  M  for  the  stack  is  analogous  to  the  one  Potter  gives  for  films  supported  by  a  substrate  and  is 
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found  from. 


where 


i; 

==  M 

I? 

Jo 

0 

M  = 


(C-4) 


is  the  intensity  transfer  matrix  for  the  beam  going  from  the  incident  medium  into  the  substrate 
and  is  a  product  of  the  matrices  in  Eqs.  (18)  and  (20)  appropriate  to  the  media, 

=  “0lCl*l2C2— 


_1_ 

3o« 


®0«  ^Oj 


(C-5) 


1 

Rfls  ^0* 

where 

(C-6) 

Tjm  comes  from  Eq.  (19),  »  3Q,/aQ,,  =  Cp^/aQ,,  Rj^  =  bQ,/aQ,  and  Dq,/ Bq^.  Similarly, 

represents  the  intensity  transfer  matrix  for  a  beam  going  fiom  the  substrate  to  the  exit  meditim 
and  is 


^0*  ”  ®».»+  lC»+  l“,  +  i,,  +  lCs  +  2"’*»  +  Nv 


(C-7) 


c.f  d.r 
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where 


^«f  “  '^»,s  +  l’^*+l.l  +  2— +  I 


(C-8) 


^sf  =  Sjf/ajf,  R,f  =  cya,f ,  and  D /a^j ,  The  intensity  matrix  Mg  through  the  substrate 

is 

Mg  =  C,.  (C-9) 

Program  STACK  uses  the  relations  in  Eqs.  (C-1  through  C-9)  to  compute  the  transmittance 
and  reflectance  of  a  light  beam  for  the  case  that  no  interference  occurs  in  the  films  or  substrate. 
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APPENDIX  D 
Complex  Arithmetic  in  TBL 

This  note  is  to  caution  users  of  Version  3.1  of  the  Interactive  Data  Language  (IDL  (tm))  for 
Windows,  from  research  Systesm,  Inc.  of  Boulder,  Colorado,  that  serious  computational  errors  may 
occur  when  one  sues  IDL  to  compute  the  square  roots  of  complex  numbers.  In  IDL,  the  real  and 
imaginary  parts  of  a  complex  number  arc  single  precision  (“FLOATING").  If  one  uses  the  SQRT 
function  to  frnd  the  square  root  of  a  complex  number  whose  imaginary  part  is  small  compared  to  its 
real  part  (or  vice-versa),  the  imaginary  (or  real)  part  of  the  square  root  may  be  in  serious  error. 

In  one  test,  IDL  program  CMPLXTST,  the  real  part  of  a  complex  number  z  is  fixed  at  one 
and  the  imaginary  part  decreases  from  one  by  successive  division  by  10.  The  program  compares  z 
with  the  supposedly  equivalent  SQRT(z*2).  Though  only  small  roundoff  errors  occur  for  imagi¬ 
nary  values  greater  than  or  equal  to  10’^,  for  imaginary  part  values  of  z  of  10”^  or  less,  the  imagi¬ 
nary  part  of  SQRT(z'‘2)  is  zero. 

To  remedy  this  simadon,  one  may  compute  the  square  root  of  any  complex  number  by  find¬ 
ing,  in  double  precision,  the  modulus  and  angle  of  its  polar  form  in  the  complex  plane  and  then  mul¬ 
tiplying  the  square  root  of  the  modulus  by  the  sum  of  the  cosine  of  half  of  the  angle  of  z  and  i  times 
the  sin  of  the  same  half  angle.  In  fact,  before  this  correction  was  applied  to  one  of  the  AEDC  pro¬ 
grams,  the  result  of  a  entain  calculation  involving  a  complex  square  root  was  1.0  instead  of  the  ex¬ 
pected  value  of  about  0.3. 
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